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Abstract: Terascale supersymmetry has the potential to provide a natural explanation 
of the dominant dark matter component of the standard ACDM cosmology. However once 
we impose the constraints on minimal supersymmetry parameters from current particle 
physics data, a satisfactory dark matter abundance is no longer prima facie natural. This 
Neutralino Tuning Problem could be a hint of nonstandard cosmology during and/or after 
the Terascale era. To quantify this possibility, we introduce an alternative cosmological 
benchmark based upon a simple model of quintessential inflation. This benchmark has no 
free parameters, so for a given supersymmetry model it allows an unambiguous prediction 
of the dark matter relic density. As a example, we scan over the parameter space of the 
CMSSM, comparing the neutralino relic density predictions with the bounds from WMAP. 
We find that the WMAP-allowed regions of the CMSSM are an order of magnitude larger 
if we use the alternative cosmological benchmark, as opposed to ACDM. Initial results from 
the CERN Large Hadron Collider will distinguish between the two allowed regions. 
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1. Introduction 

One of the most appealing features of Terascale supersymmetry JlJ is the potential to 
provide a natural explanation of the dominant dark matter component of the standard 
ACDM cosmology. This explanation is driven by particle physics motivations from data 
which are completely unrelated to the astrophysical data that motivate ACDM. In most 
supersymmetry (SUSY) models the dark matter candidate is the neutralino, although there 
are other interesting possibilities. The mass and annihilation cross section of the neutralino 
are determined by the SUSY model parameters. Adapted to a ACDM cosmological history, 
the neutralino becomes a cold thermal relic, freezing out during a radiation-dominated era 
at a temperature equal to approximately l/20th of its mass. SUSY model parameters thus 
determine the current neutralino abundance, which is usually expressed as a fraction, f2 x , 
of the critical density. 

If we assume that neutralinos constitute nearly all of the non-baryonic matter at late 
times, then their abundance is quite constrained by astrophysical data. For example, 
combining the WMAP three year data with Sloan Digital Sky Survey data on large scale 
structure, one obtains Q the (naive) 2-sigma limits: 

0.095 < n x h 2 < 0.122 . (1.1) 

Thus the neutralino relic abundance is determined with ~ 10% accuracy. 

Suppose for simplicity that we consider some subclass of Terascale SUSY models with 
the neutralino as the stable lightest supersymmetric partner (LSP). The most popular ex- 
ample is the CMSSM, inspired by minimal supergravity ||-||, where the SUSY parameter 
set is written mo, mi/2, A, tan/?, and sign(^). Imposing very rough considerations of 
naturalness, we can require mo, m^ 2 and \A\ to be less than, say, 3 TeV. If we assume 
ACDM cosmology we can scan over this parameter space, placing a mark at each point 
which would predict a relic density satisfying the "WMAP" constraints ( |1.1| ). 
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Figure 1: The WMAP-allowed points for the CMSSM, with A = 0, \i > 0, and tan /3 = 5. Points 
are shown only if the corresponding Higgs mass is > 114 GeV. The scan was performed in 2 GeV 
increments. 



The results look like Figures |T| and |2|, where we have chosen A = 0, \x > 0, tan /3 = 5 
or 30, and we only plot points which also satisfy the existing experimental bounds on the 
superpartner particles and the Higgs. The "WMAP-allowed" regions are a very small frac- 
tion of the total parameter space. More significantly, once we impose experimental bounds 
the WMAP-allowed points are not generic; for a Bino-like neutralino generic points pre- 
dict a relic density that is much too large, while for a Higgsino-like or Wino-like neutralino 
generic points predict a relic density that is much too small. Getting a WMAP-allowed 
relic density thus requires tuning to regions where various conspiracies take place, either 
enhancing the neutralino annihilation (or co-annihilation) cross section, or balancing the 
Bino-Higgsino-Wino content of the LSP. These conspiracies can be quantified by defining 
simple sensitivity measures of the predicted relic density to small variations in relevant 
SUSY parameters Q. Since conspiracy is the opposite of naturalness, such an analysis j<j]- 
Jl8| ] casts doubt upon the CMSSM, and perhaps Terascale SUSY in general, as a natural 
explanation of dark matter. 

This "Neutralino Tuning Problem" is similar to the much more famous "Little Hierar- 



chy Problem" IE, 2C] of the MSSM. While conceptually distinct, these problems are related 
by the fact that the superpartner and Higgs mass bounds from data eliminate what would 
otherwise be a generic WMAP-allowed region of the CMSSM known as the "bulk" region. 
The bulk region is generic because its only distinguishing feature is that the superpartners 
are light, as favored by naturalness applied to electroweak symmetry breaking. Thus at 
least in the CMSSM the two naturalness problems are coupled. 

In this paper we revisit the important question of whether we can quantitatively evalu- 
ate the claim that Terascale supersymmetry provides a natural explanation of dark matter. 
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This issue has not been fully resolved in the literature, due to a number of difficulties which 
will we review in the next section. We will show that a physically meaningful approach 
to this question is to compare the robustness of the connection between SUSY and dark 
matter under alternative well-motivated cosmological scenarios. 

There is ample motivation to vary our cosmological assumptions. Despite the dramatic 
success of the standard ACDM model in confronting a wide variety of astrophysical data, 
this same data is largely silent on the cosmological history of the universe before big bang 
nucleosynthesis (BBN), i.e., at temperatures above about 4 MeV. WMAP data is consistent 
with an epoch of primordial inflation, but we are very far from being able to flesh out these 
hints into a data-driven model which is sufficiently robust and detailed to make definitive 
statements about dark matter. 

A difficulty in exploring cosmological alternatives is that, because Terascale cosmology 1 
is almost a black box, it is easy to introduce new parameters which are relevant but 
unconstrained, modifying the neutralino relic density predictions by orders of magnitude. 
For example, we can dilute the neutralino density by adding entropy to the radiation 
bath some time after neutralino freeze-out but before BBN. Alternatively we can increase 
the predicted neutralino abundance by producing them non-thermally. Our solution is to 
look at alternative cosmological scenarios which are reasonably motivated, make definite 
predictions for dark matter, and are sufficiently simple and constrained that they have no 
relevant adjustable parameters. The dark matter predictions of Terascale SUSY can then 
be compared for these cosmological benchmark scenarios. 

As a first alternative cosmological benchmark, we propose the Slinky model of quintessen- 
tial inflation described in [21, 22|. This model has no adjustable parameters once we require 
that the current radiation and dark energy fractions take their WMAP-preferred values, 
and that the universe is overwhelming radiation-dominated during the time of BBN. Since 
the same inflaton is responsible for primordial inflation and for dark energy, it is not sur- 
prising that the universe is not completely radiation-dominated at the time that neutralinos 
freeze out. By making the minimal assumption that inflatons do not decay to neutralinos, 
we have a picture in which thermally produced neutralinos are diluted by a predictable 
amount, over and above the dilution of standard ACDM. This additional dilution of the 
neutralino abundance can be expressed (to an accuracy of a few percent) as a polynomial 
function of the square root of the neutralino freeze-out temperature. 

Thus for any SUSY model we can compute the predicted neutralino relic abundance 
for two contrasting cosmological benchmarks, ACDM and Slinky. The ratio of the WMAP- 
allowed regions is a physically meaningful naturalness comparison. In our CMSSM scan, we 
find that this ratio is always small. For the Slinky benchmark the WMAP-allowed points 
are more generic, without balanced mixings, mass degeneracies, or resonance-inducing mass 
relations. 

These results give a concrete measure of the Neutralino Tuning Problem in the CMSSM. 
At the same time, they define the beginnings of a straightforward program to probe Teras- 
cale cosmology at colliders. 

1 In this phrase 'Terascale' refers to the physics responsible for the WIMP dark matter, not the freeze-out 
temperature. 
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Figure 2: The WMAP -allowed points for the CMSSM, with A = 0, fi > 0, and tan/3 = 30. 
Points are shown only if the corresponding Higgs mass is > 114 GeV. The scan was performed in 10 
GeV increments, which misses a few of the highest mass points in the co-annihilation region (lower 
right). 



2. The Neutralino Tuning Problem 

This problem is most easily observed in (but not limited to) the CMSSM. We bound the 
parameter space by requiring mo, w-i/2 an d \A\ to be less than, say, 3 TeV. We remove 
regions in which the lightest neutralino is not the LSP, or in which we do not get proper 
electroweak symmetry breaking. We also remove regions which are in direct conflict with 
experiment. For most of the MSSM parameter space, the most restrictive experimental 
constraint is the direct lower bound on the lightest Higgs mass from LEP. We will conser- 
vatively write this bound as > 114 GeV. We fix the mass of the top quark to the 2005 



combined Tevatron average [23] of 172.7 GeV; our results would differ only slightly using 



the newer combined average [24] of 171.4 GeV. 

Assuming ACDM cosmology we can scan over the remaining CMSSM parameter space, 
placing a mark at each point where the computed neutralino relic density satisfies the 
WMAP constraints (|1T]). Representative results are seen by fixing A = 0, jj, > 0, and 
choosing either tan f3 = 5 (Figure [j]) , tan (3 = 30 (Figure |2|) or tan (5 = 50 (Figure |||) . The 
figures are essentially identical to those in |25|] , except that our allowed regions are smaller; 



this is because we use the new WMAP 2-sigma bounds, while 25] uses the older WMAP 
99% cl bounds, which are less constraining. 

Once we impose experimental bounds the WMAP-allowed points of the CMSSM are 
not generic. For a Bino-like neutralino generic points predict a relic density which is much 
too large, while for a Higgsino-like or Wino-like neutralino generic points predict a relic 
density which is much too small. Getting a WMAP-allowed relic density thus requires 
tuning to regions where various conspiracies take place, either enhancing the neutralino 
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Figure 3: The WMAP -allowed points for the CMSSM, with A = 0, fi > 0, and tan/3 = 50. 
Points are shown only if the corresponding Higgs mass is > 114 GeV. The scan was performed in 
10 GeV increments. 



annihilation (or co-annihilation) cross section, or balancing the Bino-Higgsino-Wino con- 
tent of the LSP. 

In Figures [l] and [2|, the lower slivers of WMAP-allowed points are made possible by 
tuning the lightest stau to be nearly degenerate with the lightest neutralino, producing 
co-annihilation at the time of freeze-out. The upper sliver in Figure is a "focus point" 



region [26], where we tune to give the LSP a significant Higgsino component. This allows 
neutralino pairs to annihilate more efficiently through t-channel chargino exchange into 
dibosons. The focus point region does not show up in Figure |l] because it produces only 
points with < 114 GeV. Apart from this, shifting tan/3 from 5 to 30 does not make a 
qualitative difference. 

For a window of very large tan (3, which is itself a tuning, we obtain results as in Figure 
|3[ Here the focus point region has expanded, the co-annihilation sliver remains a sliver, and 
a new region, the "^l-funnel", has opened up. In this region we are tuning tua — ZmLSP, 
allowing neutralinos to annihilate efficiently through an s-channel resonance. 

The Neutralino Tuning Problem of the MSSM has been known for some time, but 
this observation has not yet spawned anything like the febrile activity engendered by the 
Little Hierarchy Problem. Part of the reason is that it is not straightforward to address 
the Neutralino Tuning Problem in a physically meaningful and unambiguous fashion. This 
is most easily seen by attempting to quantify the problem. 

The neutralino relic density, assuming a ACDM cosmological history, is a physical 
observable that is not conceptually different from, say, the selectron mass. Obviously in the 
limit that the WMAP errors shrink to zero, the allowed fraction of the CMSSM parameter 
space, or any other space of relevant parameters, will also shrink to zero. The "small" size of 
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Figure 4: The cosmological history of Slinky. Shown are the relative energy density fractions in 
radiation (green/dashed), matter (blue/dot-dashed), and the noncanonical scalar (red/solid), as a 
function of the logarithm of the scale factor a(t). The early period of (exotic) matter domination 
is illustrative and has no connection to neutralinos. 



the WMAP-allowed region of the parameter space, in and of itself, therefore has no bearing 
on naturalness. It is simply a mapping from the parameter space to a particular physical 
observable with a small error bar. The same is true for more sophisticated statistical 



analyses |27||-[p9 



Another quantitative approach comes from expanding the SUSY parameter space, 
and comparing the WMAP-allowed regions for the relic density with those obtained in the 
CMSSM. But this approach has serious drawbacks. If we expand our class of SUSY models 
by varying parameters which are irrelevant to the determination of the neutralino relic 
density, we have done nothing. If we expand our class of SUSY models by varying relevant 
parameters, we are also (in general) changing the measure on "theory space", so a naive 
comparison of allowed regions is not meaningful. This makes it difficult, for example, to 
compare the CMSSM with the non-universal Higgs models defined in |3(]] , where the focus 
point and ^4-funnel regions are complicated slices through a higher dimensional parameter 
space. 

It may be possible to find a class of SUSY models which do significantly better than the 
CMSSM in an unambiguous, apples-to-apples comparison. However we would not learn 
much from this unless the new models had some theoretical or experimental motivation 
independent from our desire to explain dark matter. 

Experimental motivation could come from the Large Hadron Collider (LHC). Initial 
physics runs may provide direct access to Terascale supersymmetry. This would certainly 
give a dramatic, but not complete, reduction of the SUSY parameter space compatible 
with particle physics data. It is possible that in a large fraction of this reduced parameter 
space the predicted neutralino relic density will satisfy the WMAP (or Planck) bounds. 
While this would not answer our original question of whether generic Terascale SUSY 
gives a natural explanation of dark matter, it would make this question somewhat moot. 
In such a case emphasis will naturally shift to trying to obtain precision comparable to 
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Figure 5: The WMAP-allowed points for the CMSSM, with A 
scan was performed in 2 GeV increments. 



0, u > 0, and tan/3 = 5. The 



the WMAP errors, and ultimately to assessing whether the neutralino constitutes the only 
major component of dark matter. It will not be possible to do this using only initial results 



from the LHC. Detailed studies [31] indicate that the International Linear Collider (ILC), 
combined with the LHC data as well as data from direct and indirect dark matter searches, 
will be required to provide adequate resolution of the relevant parameters. 

In the meantime, it is prudent to assess our original assumptions about SUSY and 
about cosmology Regarding SUSY it is certainly important to consider LSPs other than 
the neutralino. However as has been already noted by many authors |2q| , without simul- 
taneously modifying our cosmological assumptions this tends to complicate the natural 
connection between SUSY and dark matter, rather than simplifying it. 

Regarding cosmology, we have made two major assumptions. The first is that neutrali- 
nos are the overwhelming dominant component of dark matter. This is not very plausible, 
given that visible matter consists of several quite different stable components. On the 
other hand this assumption also has little impact on our analysis of the CMSSM, for the 
large part of the parameter space where the LSP is nearly 100% Bino. Here we generically 
predict far too much dark matter. Thus if we instead asked that the neutralino merely 
constitute a significant fraction of dark matter, our unhappiness with the CMSSM would 
be unchanged. 

Our second major cosmological assumption was that we can naively extrapolate the 
ACDM picture back to the time of neutralino freeze-out, and that the neutralinos them- 
selves are entirely thermal relics. This assumption is plausible but not yet directly sup- 
ported by any data. The astrophysical data that supports the standard ACDM model 
is largely silent on the cosmological history of the universe before BBN. WMAP data is 
consistent with an epoch of primordial inflation, but we are very far from being able to 
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Figure 6: The WMAP -allowed points for the CMSSM, with A = 0, fi > 0, and tan j3 = 10. The 
scan was performed in 10 GeV increments. 



connect these hints to definitive statements about dark matter. 

Simply put, a top-down approach of trying to fix pre-BBN cosmology sufficiently well to 
pin down the genesis of dark matter is likely to fail in the foreseeable future. More realistic 
is to combine a top-down approach with a bottom-up approach, which uses collider data 
and dark matter detection results. Thus, at the same time that we use data to determine 
the properties of dark matter, we are simultaneously using data (some of it the same data) 
to determine the cosmology of the Terascale epoch. 



3. A new cosmological benchmark 

With the assumption that neutralinos are thermal relics, their abundance Y(X) can be 
computed as a function of X = T/m x , the ratio of the temperature of the thermal bath 
divided by the mass of the LSP. The evolution equation for Y(X) is 33]: 



§ = 5P|Sm, < „ > (Y eq (X? - YiXf) (3.1) 

where g* is an effective number of degree of freedom, M p is the Planck mass, Y eq (X) 
the thermal equilibrium abundance, and < av > is the relativistic thermally averaged 
annihilation cross section. The SUSY model determines the cross section, summing over 
the relevant annihilation and co-annihilation channels. 



We have used micrOMEGAs 2 . 0, which solves ( |3.1| ) numerically 35] and returns the 



predicted neutralino relic density as £l x h 2 . It also returns the freeze—out value of X, called 
Xp, which is roughly equal to 1/20. The program uses SUSPECT 2. 3 to compute the SUSY 
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Figure 7: The WMAP-allowed points for the CMSSM, with A = 0, (i > 0, and tan j3 = 30. The 
scan was performed in 10 GeV increments. 



spectrum from CMSSM input parameters [36|. As noted in the previous section we fix the 
top quark mass to be 172.7 GeV. 

In the Slinky model of quintessential inflation, a single noncanonical inflaton is respon- 
sible for both primordial inflation and the present day accelerating expansion. The inflaton 
potential is an exponential of a periodic function, giving repeated epochs of progressively 
weaker accelerated expansion, alternating with epochs of radiation or matter domination. 
The period is adjusted so that BBN occurs during a period of radiation domination. The 
remaining free parameter, which is the coefficient controlling the coupling of the inflaton to 
radiation (Jcq), is always small, and is turned off at late times to avoid strong phenomeno- 
logical constraints. This weak coupling is fixed by requiring that the present day ratio of 
radiation density to quintessence takes the WMAP-preferred value. The original Slinky 
allowed for the nonthermal production of dark matter via coupling to the inflaton; for the 
purposes of this paper we assume that this coupling is not present, and that dark matter 
is produced thermally entirely from neutralino freeze-out. 

Thus Slinky is a simple well-motivated cosmological model with no adjustable param- 
eters. Not surprisingly, at the time that neutralinos freeze out, the Hubble expansion rate 
H differs from what it would be in standard ACDM, since the universe is near the end of 
an epoch of accelerated expansion. Furthermore, between the time that the neutralinos 
freeze out and BBN, there is a significant entropy increase due to the inflaton coupling to 
radiation. The net effect is to dilute the neutralino relic density relative to what would be 
predicted with ACDM. Thus f2 x (Slinky) = S7 x (standard)/7, where 7 is the Slinky dilution 
factor. From a numerical analysis, we have found a heuristic formula for 7 as a function of 



7 = 0.6 + - 0MX F + O.OlXp 



(3.2) 
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Figure 8: The WMAP-allowed points for the CMSSM, with A = 0, p, > 0, and tan j3 = 50. The 
scan was performed in 10 GeV increments. 

This formula is accurate to within a few percent. The exact solution can be obtained by 
calculating the entropy production through numerically solving the evolution equations for 
the inflaton, radiation and matter, i.e. 



pe = -3H(1 + w)p e - k H(l + w)p e ; 
p r = -AHp r + k H(l + w)p e ; 
Pm = —3Hp m ; 



(3.3) 



where H is the Hubble rate, w(t) is the inflaton equation of state parameter, and pg, p r , 
p m are the energy density of the inflaton, radiation and matter, respectively. 

4. Results 



The results of the CMSSM scans are shown in Figures (g-g), for A = 0, p > 0, and various 
values of tan (3. We only show points for which the neutralino is the LSP, we obtain proper 
electroweak symmetry breaking, and we satisfy current experimental bounds including 
m h ^ H4 GeV. The yellow points predict neutralino relic densities which satisfy the 
WMAP bounds (|1.1| ). The magenta points predict neutralino relic densities which satisfy 
the WMAP bounds when combined with the Slinky dilution factor as given by (Kty. These 
points are consistent with all accelerator bounds, and we have checked using DarkSUSY3 . 14 
1 37] that the resulting spin-independent LSP-nucleon cross sections are below the CDMS 
limits H|. 

As seen in Table 1, the WMAP-allowed region is an order of magnitude larger if 
we assume the modified Slinky cosmology rather than ACDM cosmology. In the large 
tan P ~ 50 window this ratio is reduced to less than 3, but recall that this window is itself 
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tan (3 


number of allowed 
standard points 


number of allowed 
Slinky points 


5 


187 


4812 


10 


68 


435 


30 


108 


434 


50 


1423 


3578 



Table 1: Results of CMSSM scans in the mo-mi/2 plane, with A = 0, [i > 0, and various values of 
tan/3. The scans were performed in 10 GeV increments, except for the case tan/3 = 5, where the 
scan was in 2 GeV increments. The second column contains the points which give a WMAP-allowed 
neutralino relic density, as reported by micrDMEGAs 2 . 0. The third column shows the number of 
WMAP-allowed points when wc include the additional dilution factor of the Slinky benchmark 
model. 
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810 
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341 
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114.1 
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0.990 
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115.1 
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109 


492 


0.906 


1350 


170 


30 


114 


65 


1239 


113 


500 


0.887 



Table 2: Superpartner and Higgs masses in GeV for selected CMSSM points which predict the 
correct amount of neutralino dark matter using the Slinky benchmark cosmology. All points have 
A = and > 0. The particle spectra are computed using SUSPECT 2. 3. 



a (small) tuning. The large ratio for the case tan/3 = 5 is enhanced by the lack of an 
allowed focus point region for the standard case. This comparison of allowed regions is 
only weakly dependent on the size of the WMAP error bars, or on the choice of metric in 
the "theory space" . 

From a detailed comparison of the particle spectra, we find that the WMAP-allowed 
points in the Slinky case are more generic. This is not obvious from the figures, since except 
for the case of tan (3 = 50 the magenta regions are fairly close to the yellow regions where 
stau co-annihilation or Bino-Higgsino mixing is occurring. Looking at the spectra, we find 
that the lightest stau is at least 40 GeV heavier than the LSP in the lower magenta regions, 
while the Bino fraction in the upper magenta regions never dips below 0.88. The magenta 
regions do not involve the kind of balanced mixings, mass degeneracies, or resonance- 
inducing mass relations that characterize the yellow regions. Table 2 shows the particle 
spectra for some representative magenta points. 

5. Conclusion 

In this paper we have reviewed the Neutralino Tuning Problem of the CMSSM and MSSM. 
We have quantified this problem by an explicit comparison of the CMSSM for two different 
cosmological benchmark scenarios: standard ACDM and Slinky. The comparison shows 
that the standard case requires about an order of magnitude of extra tuning to successfully 
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account for dark matter. The WMAP-allowed regions for the Slinky case are still restricted 
to rather narrow strips, but this largely reflects the small (± 10%) error bar on Q x h 2 that 
we are trying to squeeze into. 

We expect that it should be possible to construct interesting new cosmological bench- 



marks from examples already in the literature p9|] -[44|. Note that a different cosmological 
benchmark with a larger dilution factor would look even more natural than Slinky. Suppose 
for example that we fix mo, A, sign(/i) and tan/?, and look at 

f{rn x / 2 ) = \dtt x /dm 1/2 \ . (5.1) 

We can regard /(m^) as a rough naturalness measure; larger /(mj / 2 ) in a region where we 
obtain WMAP-allowed points means that we will obtain a narrower strip of allowed points. 
Generically /(m^) is largest at both small and large values of mi/2, passing through zero 
at some intermediate value of m^2 where £l x is maximized. Thus a cosmological benchmark 
with a larger dilution factor will pick up WMAP-allowed points in a region where /(7711/2) 
is smaller. In addition, when we have a dilution factor 7, the relevant naturalness measure 
is not f{mi/2) but rather /(m, 1 / 2 )/7, since it is the latter quantity which determines the 
width of the WMAP-allowed region in the approximation that 7 is a slowly varying function 
of m 1/2 . 

The superpartner spectra in the standard and Slinky cases differ by at least tens of 
GeV. In addition, the Slinky spectra do not exhibit the mass degeneracies or resonance- 
inducing mass relations that occur in the standard case. Thus it should be possible to 
discriminate between these spectra using initial data runs from the LHC. This opportunity 
deserves detailed study. 

Suppose that experiments at the LHC do discover a spectrum of new particles consis- 
tent with the CMSSM, but not consistent with the WMAP upper bound as derived assum- 
ing a neutralino LSP and standard cosmology. Then we will face the interesting challenge of 
trying simultaneously to refine our knowledge of both the underlying Terascale physics and 
the underlying Terascale cosmology. Though both may be complicated, it seems prudent to 
develop a strategy that begins with simple benchmarks. The CMSSM+Slinky benchmark 
described here is a step in this direction. 
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